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Evaluation of Surface Water Resources from Machine-Processing 


of ERTS Multispectral Data 1 


I’. W, Mituiscl, W. J, Twltl, M. F, Baumgardner, R. A, Mitchell, and J. P, Cook 1 


A»r, TRACT 

Water resource data that era usoful to environmental scinntlsis 
and planners frequently arc mission, incomplete, or ohtninod ir» 
regularity, A new source of surface hydrolopical information can 
he obtained as often as every 18 days in some areas through 
mncli inn-processing of Earth Resources Technology Satellite (ERTS) 
multispectral scanner data. This research focused on the surfaco 
water resources of a largo metropolitan area, Marion County (In- 
dianapolis), Indiana, in order to assess the potential value of ERTS 
spectral analysis to water resources problems. 

The results of the research indicate that all surface water bodies 
over O.t, ha wore identified accurately from ERTS multispectral 
analysis. Five distinct classes of water wore identified and cor- 
related with parameters which included the |) degree of water silti- 
noss; li) depth of water; iii) presence of mnero and micro biotic 
forms in the water; and iv) prcscnca of various chemical concen- 
trations in the water. The machine-processing of ERTS spectral 
data used alone or In conjunction with conventional sources of 
hydrological Information can lead to the monitoring of the i) nma 
of surface water bodies; h) estimated volume of selected surface 
water bodies; Iii) differences in degree of silt and clay suspended in 
water; and iv) degree of wator eutrophication related to chemical 
concentrations, Water resources Information obtained from ERTS 
analysis will bo useful in helping to solve or bettor understand 
selected pollution, erosion, and planning problems in metropolitan 
and other environments. 

Additional Index Words: water pollution, water management, 
water quality monitoring. 


Water remurccs are vitally important to large metropolitan 
areas heea.tsc industries, commercial establishments, and 
residences depend upon the amount and quality of a re- 
gion’s water supply. Hydrological data are used to give 
insight into many aspects of water supply and quality, 
drainage and flood control, water recreation, and sewage 

'Contribution from the Department of Geography and Geology, 
Indiana Stale University at Terre Haute 47809 and the Laboratory 
for Applications of Remote Sensing, Purdue University, West 
Lafayette 1 7906, Research was sponsored in whole or in part by 
the National Aeronautics and Space Administration (NASA) under 
Grant No. Mil, 10-005-1 12, Journal Paper No, 5299 of the Agr, 
Exp, Sta,, I’nrdue University, Received B November 1973, 

^Associate Professor of Geography and Geology (ISU), Research 
Geographer (LAItS), Professor ot Agronomy and Earth Resources 
Programs Leader (I.ARS), and Graduate Research Assistants (ISU), 
respci lively. 


processing in nuisl rncinipnliLtn ureas', liowcvn, severe 
problems of djiln iicfjuisilion arc common in many lypes 
of environmental water studies. A toial or partial pain iiy 
of surface hydrological data inhibits a good regional as- 
sessment and evaluation of water resources. T his paper 
explores the use of the Earth Resources T ec Tuioiogy 
Satellite (ERTS) us a source of information which pro- 
vides regional surface water data suitable in list* inde- 
pendently or in conjunction with existing water resource 
data in order to evaluate better water quuliiy ami uses in 
metropolitan environments. It will be shown that ma- 
chine processing of I) RTS multispectral scanner data pro- 
vides a valuable method to obtain hydrological informa- 
tion that can he used In help monitor and trv; JitaU* various 
surface water characteristics such as dcplh, surlaii* area, 
volume, and turbidity .simultaneously as often as every 
eighteen days in « given study area, Theoretically over 
31 million km 1 of surface area can he scanned everyday 
by ERTS from its orbital altitude of 925 km. These 
ERTS capabilities add greatly to the ability of die en- 
vironmental scientist and regional planner to make ration- 
al water resources decisions. 

' MATERIALS AND METHODS 

Maricm County, (Indianapolis) Indiana (Fig. 1) was selected as 
an area for water resources analysis utilizing ERTS multispectral 
scanner data because: i) the county is one of the ten largest 

metropolitan areas in the United States, thus positive results gained 
from ERTS analysis will be significant to many individuals anil ii) 
the diversity of surface water features presents a good op|.(>'tumiy 
to assess the ERTS system in a hydrological application, 'Hie 
largest water hollies arc Eagle Greek and Geist Kaereoirs in die 
northwestern and northeastern sections of the county respectively. 
Geist Reservoir is the major supply of water Tor the iimaliitanis of 
tile county. There arc approximately 60 ponds or sma‘1 lakes oT 
more than 0.5 ha (1,2 acres) in the study area. Several of these 
ponds are active or inactive gravel pits, some of which are associ- 
ated with lowlands that receive overflow from adjacent creeks and 
the West Fork White River, Sewage disposal and water treatment 
facilities use water bodies (maintained by man) ailj.a cut pi die 
West Fork White River, Many ponds sir smalt lakes arc found in dc- 
pressional topographic locations r.ut associated with reservoirs, 
floodplains, or gravel pits. The width of the major river of the 
county (West Fork White River) rarely exceeds 150 tn, while the 
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average width of smaller tributary ritrcumtt (be. Fall Greek) is less 
ill, m (id in. 

Frau haitils of hightpuliiy digitized, multi spectral data were ob- 
tained firm .111 I, ll i.S |i,is% over (In' study arr.i on .'ll) Scplcmtivr 
l')?» (It.mn •), , r i ,(i /in; |t, mil fit 7 pm; |l,iml II, ,?-,K pm; anil 
It, mil 7 , ,N I . i pm). Computer implemented priH fxsIilg of spectral 
ilii.i was used lo assess ilu* iliKiiilmliuti .mil tlnir, mi'll, ilirt. of stir- 
f.uv water in Marion County. .Selected implication* of (he research 
results obtained from Ilu* liRT.S hydrologic analysi.s ,iru explored. 

The first step in the analysis was lo generate « gray scale line* 
printer map of the couniy in •hunrl C (0.7-0, 8 pm), Histograms 
dial Imd been calcnl.iled |irevlmo>ly allowed die map lo show six- 
teen data ratines, each of which was represented iiy a different 
alphanumeric symbol. Hand 0 was selected for display because of 
the high probability dial die lowest spectral reflectance values 
{darkest tone) in dial portion of the spectrum would lie water, 
Ilu* printout was compared with color air photography ({round In- 
formation (Indianapolis Power and I-iph i Company, Indianapolis/ 
Marion Comity Color Airphoto Mosaic, 1072) to determine which 
of those dark tone areas wc*” water and which were cloud shadow. 
Aftcrlocating the water be ' die county on the gray scale map, 
line and column eoordin-.es were chosen and recorded for rec- 
tangular training sample areas of water, because of die irregular 
shape of the reservoirs, u number of samples were selected, 1G from 
Gelsl Reservoir and seven from Ragle Creek Reservoir, A single 
sample was also chosen from each of ten smaller surface Water 
babies. Care was taken in selection of die 33 samples In order to 
obtain a representative sampling from all of the area’s surface water. 
The spectral responses of these water samples in each one of the 
Tour ERTS bands were used to help train a computer to Identify 
automatically water everywhere throughout the study area. 

These samples were submitted to a clustering algorithm program 
(U'atkcr, A. G,, and D, A. l.andgrebe, 1071. A minimum of dis- 
tance approach to classification, I.ARS Information Note 100771, 
p, 129-133. Purdue University) which requested delineation of the 
spectral responses of water data Into five spectral classes Using all 
four bands of data. The output from the clustering program was 
maps of each of the 33 sample areas which Indicated the duster 
class of all points (each point Is an area of approximately GO m by 
30 m and is known as a Remote Sensing Unitor RSU) in the sam- 
ples by different alphanumeric symbols, 'flic characteristics of 
four of the clusters (each cluster represented a distinct spectral 
ciass of water) indicated that they had Informational value, hut the 
fifth was considered to be a non-informationul (non-water) spectral 
class. Training samples representative of each one of the four 
water classes were delineated directly from the cluster maps, and 
were used ultimately to train a computer to identify automatically 
different classes of water. 

It was not advisable at this stage in the analysis to classify all of 
the data points in the county, because training samples were not 
yet defined for other types of land cover. Some of these features, 
particularly cloud shadow, could easily be miseiassified ns water. 
Thus, training samples were chosen for the other land cover types 
In the county, including single-family (newer) residential, multi- 
family (older) residential, grassy (open) areas, trees, commercial/ 
Industrial, cloud, and cloud shadow. Statistics of the relative 
spectral responses (means, standard deviations, and covariance 
matrices) were calculated for the four water classes and the seven 
other classes. All of the data points in the county were then clas- 
sified by a Gaussian maximum likelihood classifier (Phillips, T. H. 
1973. i.ARSYS users manual. LARS, Purdue University) into one 
of the eleven classes. The results were displayed by a lineprinter 
using different alphanumeric symbols for each class. Figure 2 con- 
tains a synthesized photographic display of the complete (1 ! class) 
classification of Marion County, 

RESULTS AND DISCUSSION 

Virtually 10(1% of all standing; water bodies over 0,5 
ha were identified' in (lie HUTS classification analysis. 
The .shape and surface area of the large water bodies 
(particularly the two reservoirs and to a lesser degree the 
largest ponds) were determined very accurately from the 
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Fig. 2— Digital images of Marion County, Indiana land use clavufit., 
tion. Each one of eight (seven non-wotor and one water) t.lav.i", 
of phenomena identified In Marion County through machine 
processing of ERTS spectral data was displayed using an .iv>igrnti 
gray-scale tone on a cathode ray tube (CRT). All cluv.es of 
water were combined Into one class to faeditatn photographic 
clarity; however, water subclasses arc possible to obtain to th.s 
image form. Figure 2a encompasses all of Mar mo County ami 
small portions of adjacent counties (!> 900 km J ). The t.irgest 
water bodies, Eagle Creek and Gcist Reserve rs. are s rr ■ « 
In the western and northeastern portions of ih» image. Fi-,,ir. , 
2b and 2c are digital display enlargements (A X) of the. Deist ,md 
Eagle Creek reservoir areas respectively. Figure 2d »s a maxi 
mum site digital display enlargement (16X) of the Gtv,t 
Reservoir area. 
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IK IS i.Mijuiiii in, pit nicnii'i! 1 1 , i-isilii.il imt ill [lit* hltiilv 
n> i lln Ii.ipi uni mii l.n r . 11 * .1 <il I hi' smaller Mantling 
" ni l luMti v Ui ir less .i< c ii i,i li ly (Irii imiiicil, .il t hour'll 
me nli'nMii ii i ■ >n ul ihr w.iln lii«rli’*s themselves w is 
in . nh I * iii",' ,h * in-, lie, Uisiuiiiiiii ill shape ,iml area nf 
« -it*’* In'll,' i iifi.iilv in* leases .is iht' si/e of llu* water 
"•i ' I* • ,i ; • ,i* .l iliM i eases. 

Kim h is Mime ill. in (ill in wide Were generally iili'ilti* 
sits! in I K i-S dialysis, Sin'. mis nr rivers 30 to (it) m witltr 
Mi tt' ii!* unit, tl in t .isimi.illy, while those under MO in wide 
u't ie imt illumin'd. Consequently, due in the resolution 
1111111’ tit I.R I S d.i 1 . 1 , only the West Fork While River wits 
< (insistently itlentiliitl, althmij'h portions of smaller run* 
ninq svali'i features were classified as water. 

l out spi'i t rally distinct classes of water were identified 
in M.inmi Comity (water ! through water 4). A fifth 
* lass ol water (referred to as water 5 throughout this dis. 
iiissiim) was not spectrally identified, hut nevertheless 
seems in esist ,is part of the class rotntnmf. Variations 
in the parameters that aliened the spectral response of 
was* i pheiioincn.i have resulted in a paliern of water class 
disiiiimtion that is not random, The distribution of the 
live classes' nl water helps to give insight into ihe param- 
eters that alleel the spectral responses of water features 
that wete ohi, lined limn the four hands of RUTS data. 

The areal distribution of the spectral classes of water in 
Marion County was complex. Water 1 through water 3 
were found in a similar sequence in both tingle Creek and 
Ceist Reservoirs, The extreme northern section of each 
reservoir was characterized by water 1; the north-central 
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Fig. 3-Alphanumeric representation of Gcist Reservoir water 
classes a-, determined from computer-implemented classification 
of water using four EP.TS spectral bands. Each alphanumeric 
symbol is equivalent to an area of approximately 60 by 80 m. 


net lion contained watvr 2; and mosl of the remaining 
area was water 3 (figs. 3-4). Wafer 5 was oflin lonrid 
adjacent to the reservoir shoreline. Virtually no waler4 
was found in reservoir areas, 

All four spectral classes (and the partial water thns 
water 5) i hnrat Ifiiml the waler assot inlctl with Mantling 
waier bodies (excluding reservoirs), Frcrpu nily .landing 
water hotlies were dominated by one class of wafer only* 
while in other cases two nr three t lasses ol wafer were 
more characteristic. Approximately 15%, 25%, 20%, and 
40% of the standing water bodies (excluding reservoirs) 
were dominated by water 1, water 2, water 3, and water 
4, respectively (Fig, 5). The distribution of the lakes and 
ponds dominated by a given water class appeared to he 
nearly random. 

Rivers that were sufficiently wide to permit the KR’1.3 
identification of a spectral class of water were always 
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Fig. 4— Alphanumeric representation of Eagle Crock Reservoir 
water ctassis as determined from computer-implnmonlcd classi- 
fication of water using four ERTS spectral bands. Each al- 
phanumeric symbol is equivalent to an area of approximately 
60 by 80 m. 
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identified .is waiet I (Fig. ft). Very narrow (under (it) m 
miming water bodies wcic I'rcf iiichi ly identified as water 

| In* Kriiion- Sriuiiitii Umis m USD's lli.il i, in bo urn* 
•tili'i i'll w itci r. Ii.im* .1 |i,iitnii nl meal distribution lli.it is 
ili.inu u’iisiii ol w.itri .uni they me I omul in zones when’ 
w.un .nut iiuii w.iU't iiiU'i l.in’. This type nl Witter is 
il.issilieil us i tunmene {tine til the seven mm-water 
classes), Imt it tun be distinguished accurately I'mm actual 
commerce in the study men through u spatial analysis nl" 
speur.il data ih.it is combined with uvnitniile ground ob- 
servation infunmilion. For example, .1 meandering I-J’.bu 
wide linear put tern of commerce in a rural area and a erne 
RSU layer of commerce adjacent to an isolated part of a 
reservoir are incongruous, hence these patterns can he 
identified as Water 5 .uher than commerce. There are a 
lew areas (primarily within the more highly-urbanized 
sections of the country) where the spatial and spectral 
characteristics of water 5 are more difficult hut possible 
to differentiate from commerce. 

Spatial analysis of the spectral class commerce identifies 
I -HSU wide hands oT water 5 along several narrow (< 00 
m) watercourses, adjacent to the hunks of wider (> (if) m) 
watercourses, and on shorelines o, .Hiding water bodies. 
Apparently, water fi is a mixture of water reflectance com- 
bined with relict lance from adjacent non-water phenom- 
ena (i.e,, soil, vegetation, roads, and other features, found 
on both sides of a narrow river or simihr features along a 
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Fig. 5-Alphanunncnc representation of water ctasses of smalt lakes 
and narrow stream (the stream is indicated by the winding line 
of symbol Y m the wester's parr of the figure) as determined 
from computer-implemented classification of water using four 
ERTS spectral bands, Each alphanumeric symbol is equivalent 
to an area of GO by 80 m. 


sIlori’lilK’ in lilt* cast’ of a standing wait r h< »>f > |. Ii may in 
feasible to separate accurately a (lass wain a frmn/nw 
merer solely through spectra! analysts |,«.v.rvi r, ai i!u. 
lime the use of spatial patterns nl tin >.f .* > ii.d < I t . , mu 
merer supported by availafde gromm ,ii*s« tvaiuni 
lion provitles the most act orate method ol dc Ion .iting a 
fifth class of water in the study area. 

Il is important to attempt to identify w m i r < .is a -.up 
plcmenl to inform. ilion gained from analysis ■>! waur t 
through water 4 heeanse: i) improved at t m.n y in t -.tiroa 
lion of the surface area of water bodies is attained: uj 
more complete and accurate idenlificalion ol very small 
water bodies (< 1.0 ha) isassured; iii) many rivers iii.it an 
of insufficient width to he spectrally identified as w.u. r 
arc subject to identification as the partial water i lass 
water 5; anti iv) .semipermanent cpiasiswarnpy aims t\ >'i 
and vegetation forms that periodically arc iulcrlanri with 
shallow standing water (luring periods ol high uiof.iurcq 
that generally arc impossible to identify as a spe< n.tl t la .s 
of water, arc subject to water 5 identification. Water .1 
cannot be used to analyze biological, chemical, or physii d 
characteristics of the water because too many strong non- 
water influences affect the spectral reflet lam e ol i-.uh 
RSU. 

Ground observation information, upon width hypothe- 
tic* in this study were based, was obtained lor the study 
area water bodies from maps, (Indianapolis and Marion 
County Map, 1072. Cram Company, Indianapolis}, acred 
photographs, published reports, (Division o! Planning and 
Zoning Staff. 1072-73, Water quality i ontml program 
Work papers 3-G, 9-12, and Water quality management 
plan— Indianapolis Metropolitan Area, Department ol 
Metropolitan Development, IndhinapoIivMarioti County, 
Indiana), personal interviews, (R. M. Holding, Senior 
Planner, Indianapolis Department of Metropolitan Devel- 
opment and D. Ifochcm, Engineer, Indianapolis Wain 
Company, Personal Communication), and visual ub-aivu 
lions, Equal amounts of data for each water hotly were 
not possible to obtain. Sonic potentially impoii.mt 
ground observation data could not be obtained; however, 
the sum total of all data collected permitted an evaluation 
oT several parameters that affected the spectral churn ici- 
islics of water. The information gained from this analysis 
will be. valuable to establish ground observation proto - 
dttres for future Marion County ERTS liydiologu.il 
studies. 

The parameters that appeared to have the greatest in 
fluence on the spectral characteristics of water are: i) 
variability in water depth; ii) degree of water turbidity 
particularly related to the amount of silt or clay particles 
in suspension; iii) distribution and amount of m.urobn 
tanical forms (aquatic weeds); iv) distribution and amount 
of micrnhotunicul forms (algae blooms); v) chetnii .il im 
purities; and vi) influence of adjacent earth surface lea 
lures that are not an intimate part id a water body (i.e., 
tree adjacent to a shoreline). Each different class ul w.iiei 
was influenced by varying combinations of these major 
parameters, The results of the water t lassification an iKsis 
arc presented in summary form ( fable 1 

An analysis of the data presented in fable 1 indicates 
that water 3 and water 4 were the classes least modi tied 
spectrally by non-water phenomena, while the spec tr.d 
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( h.ir.u icristii s ill' water ft reflected most strongly non- 
w.iut U’.mni s. Water i »f ilw* »!>»•« ti .il water < lasses was 
most noticeably affected spectrally by non water niateri* 
ais funnel in i lie water. High turbidity causa! by mixing 
of silt and clay (via running water, wave action, 01 human 
disturbance ), particularly in shallow water, was the most 
common condition associated with water 1. High biomass 
content, in the water (with very little silling) characterize 
some areas with water 4. Ponds and lakes that were in- 
fluenced strongly by human activity frequently are prone 
to a water I classification. Many water bodies with water 
4 were in ipiciils of inorganic and organic materials 
(sewage, lertilizer, and some sill) obtained from nearby 
residential areas and farms. It is speculated by the authors 
ili.il water 1 is often more nutrient rich overall than are 
the other standing water body classes. A mi 1 1 lent rich 
status in the water rowhf have induced above average 
biotic activity m these areas. Several gravel pit lakes also 
havcawaiei 1 classification. In these areas a combination 
of gir.n depth and moderate biotic activity (with little 
silting) affects the spectral characteristics of the water. 
Water 2 appeared to be transitional between the water 
classes strongly influenced spectrally by water containing 
suspended solids and water that was relatively lrcc of 
solids. 

The ohseivatioiis made between water classes and water 
characteristics 'Table 1) are supported by electromagnetic 
wavelength reflectance theory. The water class with the 
least silt and other solids should have the lowest spectral 
reflectance because of the great spectral absorption in ail 
ERTS hands. Water mixed with aquatic green vegetation 
characteristically should have low spectral response also 
(albeit higher overall than pure water) because of the 
ability ol green vegetation to absorb light moderately to 
strongly in tin- visible wavelengths and for water to absorb 
light very strongly in all ERTS bands. Water that con- 
tains large amounts of silt and clay (and possibly chemical 


precipitates) has relatively high rclieri.mre in the visible 
part of the spectrum because of the high spcdr.il reflet- 
tiviiy of the sill and clay fraction in the water. 

TJieorelknlly, water 1 (high amount of silt and other 
solids possibly present) should have the login st overall 
spectral response of the spectral wider classes. Water 2 
(moderately silty with biotic influences possible also) 
theoretically should be lower in overall sprt tr.il response 
than water 1, but higher in spectral response than water 
II (less silty than water 2 with biotic influence possible). 
The lowest spectral response overall of the spectral water 
■ lasses theoretically should be water ! which frequently 
contains aquatic vegetation (low reflectance) and very 
little silt and day in suspension. The actual relative spec- 
toil water classes generally agree with these theoretical 
expectations. In theory, the highest spectral responses of 
all classes that arc influenced by water is water 5. The 
spectral responses of water 5 arc determined more by sur- 
rounding relatively highly-refler live non -water phenomena 
titan water itself, thus its overall spectral response is high- 
er than that of any of the lour spectral (lasses of water. 

These explanations of spectral variations in water i lasses 
undoubtedly are valid for a very large majority of all 
water bodies of Marion County. A dearth of data for 
some of the known parameters that cause spectral varia- 
tions in water coupled with the possible cxi fusion of 
other parameters that were not considered it. tins analysis 
results in conclusions that arc subject to slight modifica- 
tion, An analysis of the available data does provide work- 
ing hypotheses which can be tested, modified, and im- 
proved as future ERTS data are analyzed for their by- 
drologieal content. 

CONCLUSIONS 

It has been shown that computer-implemented pro- 
cessing of ER'i’S data lias a potential use for hydrologit at 
studies of a metropolitan county. AU voiding water 
bodies greater than '0.5 ha are relent .ily, thereby 

permitting an accurate estimation of tie * t il .ticanl sur- 
face water in the county. These area data when com- 
bined with topographical information will permit an esti- 
mation of volume as well us area of the surface water re- 
sources of the country. Temporal variations in area and 
volume of surface water are available theoretically every 
18 days. 

A good possibility exists to differentiate lie tween water 
relatively free of suspended sill and day from water con- 
taining a large silt and clay load. Monitoring this form of 
turbidity in water bodies can be useful in accessing sources 
and intensity of erosion in the study area, lit selected 
cases, this information could be used to help make region- 
al planning decisions. 

The degree of water eutrophication, as measured by the 
vegetation content of water bodies, may be possible to 
moniior using ERTS technology. These data arc im- 
portant for water pollution studies that attempt to trac e 
and analyze the effects of nutrient -rich effluent on the 
environment. 

Planning water use for recreation and residential pur- 
poses requires water quality data. The cxuri pragmatic 
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natine nl I.UIS sprit ral analysis applied In wain pmb* 
lems iviu.iins in In* ilnmmsiialnl, .mil in pari depends im 
furtlni MM'S hydrn|ngi< al analysis nl a si inly arm al dif- 
IVieni utm'Mil i hi* ymi, However, irsnlis I mm this initial 
analysis an* veiy i'lii (imaging because Ihey indicate tltal 
significant pi.nmal applications nl spet tral data arc fetish 
hie lu nsc in help evaluate water resources. 

Ait IK-d.iy cycle monitoring the water resources is un- 
likely because nl changing weather rundilions, lull, fre- 


tpicnl ev.diia.'inir nl w.iicr re-*- nines tlt.o would la vain 
aide Lu niuiny planning and i nvimmin dial nllm.if c. 
ihle. It is likely that irinniinritij', nl « uitipuii r imp!* 
,'d processing of MMS spnlr.il data will la < m* *.l 
trust effective melltrnls m obtain a li t li d type. ol 
, .-<|ii.ility small-scale hyrlrnloijie .il data. I in' ttf i im.iit- 
elfectiveness uf monitoring will tlepemf on tin- afolity in 
understand heller how various parameters .dint tin- 
spectral responses of water in the study area. 


Column Studies of Soil Clogging in a Slowly Permeable Soil as a Function of Effluent Quality 


1 


T. C. Daniel and J. Domna 1 


ABSTRACT 

Clogging ns a function of effluent quality was Investigated In 
cores of the very slowly permeable Almcnn silt loam soil which 
offers problems for conventional on-silo liquid waste disposal. Un- 
disturbed GO cm long cores were subjected for approximately 120 
days to constant ponding with simulated septic tank effluent, ex- 
tended aeration effluent and distilled water. Column Influents and 
effluents wore monitored with respect to chemical oxygon demand 
(COD), biochemical oxygen demand (BOO), and solid residue frac- 
tions. Column influents differed markedly In COD and BOD con- 
tent hut column effluents had consistently low contents indicating 
the high runovativo capacity of the sail, In situ tensiometric, redox, 
and flow rate measurements indicated development of the most 
severe barriers to flow in columns ponded with tow BOD aorntod 
effluent, followed closely by those ponded with high BOD septic 
tank effluent. No Hjrriors developed in columns ponded with 
water. Total concon. rations of solid residue fractions In tiro two 
effluents and the cumulative load of solids applied to the columns 
did not differ significantly, but particle sizes in the aerated effluent 
were smaller. Increased pore clogging in aerated Influent treat- 
ments points to the significant role of effluent solids in the clogging 
process in slowly permeable clayey soils. Additional studies are in 
progress to better define critical waste characteristics as related to 
soil clogging. 

Additional Index Words: septic tank offluont, extended aera- 
tion, soil disposal, liquid waste, 


Slowly permeable soils, stub us the Almcna silt loam 
(Aerie Glossuqualf) unci other soils with comparable per- 
meability characteristics comprise approximately 800,000 
hectares (2 million acres) in Wisconsin and are defined as 
unsuitable by the current Health Code for on-site subsur- 
face disposal of septic tank effluent due to the low per- 
meability. Measured percolation rates of 27.5 inin/cm 
(70 min/huh) in die topsoil and 39.4 min/em (100 min/ 
inch) in the subsoil exceed the critical rate of 23. (i min/ 

'Contribution from the Soil Science Department, College of 
Agricultural and l.ifc .Sciences, Univ. of Wisconsin ami (lie Gcol. 
Nat. liist. Survey, Univ. extension, Madison, Wisconsin. This re- 
search was part of die Small Scale Wjsic Management Project, Univ. 
of VVistot.jin funded by die Upper Great Cakes Regional Commis- 
sion and die State of Wisconsin. Received 9 Dec. 1973 

2 Assistant Professor and Associate Professor of Soil Science, re- 
spectively. 


cm ((if) min/inch) (24, 25). The rapacity of a soil to » 
cept and conduct liquid can he heller expressed by mu 
sidcring hydraulic conductivity data width is phy ,i* .dly 
well defined (7), Measurements in si l ti in the Ainu im silt 
loam yielded Af Si , t values of 4 an f thy (0,98 gnl/fi 2 per 
day) for the topsoil and 2 cm/day (0,49 gal/ft 2 pn d.tyj 
for the subsoil. 

The low hydraulic conductivity of Almetia silt loam v.ii 
can he illustrated by considering the derived themt in ,d 
loading rate as a function of the percolation rate at (Hid- 
ing to O = 5 // 1 (Q loading rate in gal/lt 1 per day; t 
percolation rate in min/inrii) which is the standard equa- 
tion. used to size seepage beds (24). A percolation rate nl 
39.4 min/cm (100 min/inch) translates into a loading rate 
of 0.5 gal/ft 1 per day, or 2.04 cm/day, which is equal to 
A' S ai, the conductivity at saturation, of the horizon, 

Crust formation and clogged layers within diuitijgi- 
beds will appreciably reduce soil permeability (7). Clog 
ging of soil pores by suspended particles from diluent 
may decrease infiltration (11). Other data suggest tli.it 
clogging may he due to production of gums, derived Irmn 
organics in the liquid waste and m the soil pores 11. 

14, I5( 16, 17,23). 

Prior research iias often assumed ih,u < logging is id. on- 
ly related to the carbon load of the diluent, It SmsIhmi 
• suggested, therefore, that aerobic treatment oi waste 
water which can substantially reduce the i.ub*ni load 
could lengthen the life of the system by rnlm an; the 
clogging problem ( 1 2). Investigations have then Inn into 
limited to determinations of ambient l.u tot-- wltu ii nnhu e 
or alleviate the clogging problem within the seep., ;, i.nl. 
Unfortunately, research has lu-en to i ruled to inlmun 
studies using artificially aggregated soil till materials f 1 « i . 
rather than undisturbed cores in which flow patterns. >.l 
liquid arc significantly different (2). in this studs large 
undisturbed cores of the Almcrta silt In. on w» re exitai ;rd 
in the field and pretreated stall ti nt minimal tiisimb.mi 
of structure and soil solution occurred pn.i'r r»- , • . -e . 
column studies, The formation ol flow inliiinitng l.«yns 
as a function of carbon load was billowed with tensiom- 
eters. 
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